Drop impact upon micro- and nanostructured superhydrophobic surfaces 
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We experimentally investigate drop impact dynamics onto different superhydrophobic surfaces, 
consisting of regular polymeric micropatterns and rough carbon nanofibers, with similar static con- 
tact angles. The main control parameters are the Weber number We and the roughness of the 
surface. At small We, i.e. small impact velocity, the impact evolutions are similar for both types of 
substrates, exhibiting Fakir state, complete bouncing, partial rebouncing, trapping of an air bubble, 
jetting, and sticky vibrating water balls. At large We, splashing impacts emerge forming several 
satellite droplets, which are more pronounced for the multiscale rough carbon nanofiber jungles. The 
results imply that the multiscale surface roughness at nanoscale plays a minor role in the impact 
events for small We < 120 but an important one for large We > 120. Finally, we find the effect of 
ambient air pressure to be negligible in the explored parameter regime We < 150. 



INTRODUCTION 

Drop impact on a solid surface is ubiquitous and cru- 
cial in a variety of industrial processes. The interplay of 
several experimental conditions often makes the impact 
dynamics surprising and too complex to elucidate [l|. 
The control parameters include the impact velocity, liq- 
uid density, viscosity, surface tension, and the wettability 
and roughness of the substrates. In the past decade, the 
investigations have actively focused on the impacts upon 
superhydrophobic substrates [3, 0, 0, H, @, 0, @, due to 
their emergence and opportunities in a wide range of ap- 
plications: for instance, self-cleaning and anti-coating for 
lab-on-chip devices, ink-jet printing, spraying techniques, 
and coating processes. The most commonly studied sub- 
strates are silane-coated microstructures and hydropho- 
bic microtextures. These samples, possessing a combi- 
nation of chemical hydrophobicity and physical ordering 
roughness, display superhydrophobicity-with large static 
contact angle above 150° and small contact angle hystere- 
sis within 5°. 

In this study, we investigate the dynamics of drop im- 
pact not only on superhydrophobic microstructures of 
controlled roughness but also on multiscale rough sur- 
faces of carbon nanofiber jungles (CNFJs). The aim is 
to compare the impact behaviors upon both types of su- 
perhydrophobic surfaces and to examine the effect of the 
details of the roughness. Moreover, complementary to re- 
cent studies of drop impact using superhydrophobic sur- 
faces, our work examines the pressure of the surrounding 
air, as an additional control parameter, which recently 
has been shown to affect drop impact upon dry wetting 
surfaces 0, Ell, E3 • 

Chemical and/or geometrical surface modification af- 
fects the solid-liquid interactions. This offers a conve- 
nient control in surface physics. For example, physi- 
cal roughness increases (decreases) the contact angle for 



hydrophobic (hydrophilic) substrates. Recently, super- 
hydrophobic surfaces through controlled microstructures 
(see Fig. QJ) have been developed, advanced, and in- 
vestigated for their water-repelli ng H , Hjj . Q3] and hy- 
drodynamic slip applications [H, [l6f 7 A water droplet 
when deposited on a micropatterned hydrophobic sur- 
face can exhibit a metastable state of heterogeneously 
wetting with air trapped between the liquid and surface, 
forming a "Fakir" drop [17|, LL8[ . A spontaneous transi- 
tion can occur from "Fakir" to homogeneously wetting 
"Wenzel" state with a smaller contact angle, even in the 
case of zero impact velocity [19, 20]. At the breakdown of 



superhydrophobicity, intriguing filling dynamics of water 
infiltration was found to possess multiple timescales and 
to tune the shape of the wetted areas by the geometric 
dimensions of the micropillars [2l|. In drop impact ex- 
periment, the threshold of the impact velocity V c mark- 
ing the transition from Fakir to Wenzel state is found to 
depend on the geometric dimensions of the hydrophobic 
microstructures 0, 0]. So far, geometric arrangements 
and the impact velocity have been the main control pa- 
rameters for the drop impact onto superhydrophobic sur- 
faces [H, Q . The change in surrounding air pressure has 
not yet been explored in this context. We will also study 
this effect here. 

The analyzed superhydrophobic substrates consist of 
carbon nanofiber jungles. Carbon filaments are formed 
catalytically in metallic catalysts, particularly in Ni, Fe 
and Co based catalysts, used for the conversion of carbon 
containing gases, e.g. in steam reforming of hydrocar- 
bons and Fischer- Tropsch synthesis [22[. In these cases, 
the carbon filament formation was detrimental for oper- 
ation as they plugged reactors and deactivated catalysts. 
Fiber type carbon nanomaterials can be classified into 
three types, namely, Carbon Nanofiber (CNF), Carbon 
Nanotube (CNT) and Single Walled Nanotube (SWNT). 
In CNF the graphitic planes are oriented at an angle to 
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the central axis, thus exposing graphite edge planes. In 
CNT the graphitic planes run parallel to the central axis, 
in this state only basal planes are exposed. CNT is also 
referred to as Multi- Walled Carbon nanotubes (MWNT) 
or parallel carbon nanofibers. If the fiber consists of only 
one graphene sheet that is oriented in the direction par- 
allel to the fiber axis, it is called a single- walled carbon 
nanotube. We refer the reader to Ref. [23|, Hil, HH for the 
history, studies, and reviews of catalytically generated 
CNFs. 

The diameters of CNFs, in general, range from a few 
to hundred nanometers and their lengths vary from mi- 
crometers to millimeters [23[. The novel physical and 
chemical properties of CNFs include high surface areas, 
large elasticity, and low electrical resistivity and hence 
make CNFs good candidates for catalysts, catalyst sup- 
ports, and selective adsorption sites [25[. Particularly, 
the highly porous layers of entangled nanofibers (CNF 
jungles) are suitable as catalyst supports for liquid phase 
reaction, allowing fast mass transfer to catalytic active 
particles deposited on the CNFs 0, [2?]]. In this work 
entangled layers of CNFs are grown directly on iron and 
stainless steel metal foils, and we experimentally study 
the dynamical effect of water drop impinging upon these 
CNFs. 

EXPERIMENTAL SECTION 
Superhydrophobic samples 

Our superhydrophobic substrates are catalytically syn- 
thesized CNFJs and micropatterned polymers. Fig. [T] 
shows representative scanning electron microscope 
(SEM) images of used samples. As revealed by Fig. [T] 
(a) and (b), catalytically generated CNFJs are composed 
of entangled and bundled carbon nano-filaments forming 
uncontrolled, microscopic roughness. In contrast, Fig. [T] 
(c) shows the top-view of an ordered, microstructured 
substrate comprising micro-pillars periodically placed in 
a structured lattice. 

Synthesis of Carbon Nanofiber surfaces. The 
carbon nanofibers (CNFs) were produced by catalytic va- 
por deposition (CVD) from carbon containing gases using 
a metallic catalyst [25| . Two types of catalytic materials 
were employed: they are iron (99.99%, Alfa Aesar) and 
stainless steel Type 304 (Fe:Cr:Ni 70:19:11%, Alfa Aesar) 
foils of 0.1 mm thick. Round samples of these metal foils 
(10 mm in diameter) were prepared by an electric dis- 
charge wire-cutting machine (Agiecut Challenge 2, GF 
AgieCharmilles). The foils were degreased ultrasonically 
in acetone and dried at room temperature before loaded 
into the CVD chamber. In addition, hydrogen and ni- 
trogen (99.999% purity, Praxair), and ethylene (C2H4, 
99.95% purity, Praxair) were used for the CNF forma- 
tion without further purification. 




FIG. 1: (a) and (b): SEM pictures of two representative used 
samples of carbon nanofiber jungles, showing uncontrolled, 
multiple-scale roughness. The porous fiber-bundles span a 
few tens microns and the fine fibers are submicron thick, (c) 
A SEM top-view image of a representative microstructured 
substrate consisting of 5 /am wide pillars, yielding a well- 
defined roughness by the geometric regularity [28j] . 



The CVD reactor consists of a vertical quartz reactor, 
with a porous quartz plate centrally placed to support the 
metal foils. The temperature was raised from room tem- 
perature to 600° C at a rate of 5K/min. The samples were 
first pre-treated in a hydrogen/nitrogen mixture (20% H2 
with 80% N 2 ) with a total flow of 100 ml/min for 1 hr 
at 600° C. After the pre-treatment, ethylene was fed into 
the reactor (20% with N 2 /H 2 , 20%/60%) at 600°C for 
2 hrs (with stainless steel foil) or 3 hrs (with iron foil). 
During the heating process at the high temperature, the 
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decomposition of C2H4 led to the growth of carbon fila- 
ments upon the metallic catalysts. The concentration of 
hydrogen and the total flow were kept the same during 
pre-treatment and deposition. Finally the ethylene and 
hydrogen gas streams were shut off and the whole system 
was cooled down to room temperature under nitrogen at 
a rate of 10 K/min. 

Preparation of the polymeric microstructures. 
Micro-patterned substrates are composed of PDMS elas- 
tomer (Polydimethylsiloxance, RTV 615 rubber compo- 
nent A and curing agent B, GE Bayer Silicones). The 
fabrication of precise and controllable microstructures 
was achieved via a micro- molding method [2^]. This 
technique is to cast a polymeric film from solution on a 
molding wafer of desirable micropatterns. The PDMS 
films were obtained by mixing the rubber component 
A with the curing agent B (10 : 1 wt/wt). The mix- 
ture was degassed and then poured onto the mold (~ 1 
mm) and cured in an oven for 3 hrs at 85° C. The cast- 
ing molds consisted of diverse arrays of micro-patterns, 
which can produce periodically arranged, round or rect- 
angular polymeric micro-pillars of height h, diameter or 
width w, and the interspacing a. With different lattice 
arrangements, ft, a, and w, we controlled the roughness 
of the polymeric microstructures. We characterize the 
surface roughness 1Zf by the ratio of the total surface 
area to that projected on the horizontal plane For the 
experiments presented here, w = 5 /im, ft = 6, 10, and 
20 /xm, and a is varied between 1.27 /im and 5 /im. The 
corresponding TZf spans between 1.5 and 9. 

Contact Angle Measurements. The contact an- 
gle measurements were performed with a milli-Q water 
droplet of 4 /il, matching the volume of impacting droplet 
of ~ 1 mm wide in radius, with the Laplace- Young fitting 
method (OCA 20, DataPhysics). Catalytically generated 
CNFJs exhibit random topography of multiscale rough- 
ness, and thus we need to ensure the homogeneous wet- 
tability of such surfaces when performing drop impact 
experiments. The static contact angle Osca of CNFJs 
reported here was measured at different locations on the 
same sample, both before and after the drop impact ex- 
periment. The scattering of these data is used to estimate 
the error in Osca- 

Drop impact experiments 

The general experimental procedure consists in releas- 
ing an impinging milli-Q water droplet from a fine needle 
(0.1 mm inner diameter) with a syringe pump (PHD 2000 
Infusion, Harvard Apparatus) at different heights to vary 
the impact speed upon the solid surface. The balance be- 
tween the surface tension and droplet gravitational force 
sets the droplet size, which is « 1 mm in radius within 5% 
deviation. The whole experimental setup was enclosed 
in a chamber, connected with a vacuum pump so as to 



control the pressure of the surrounding air. The dynam- 
ics of drop impact was recorded by a high-speed camera 
(Fastcam SA1, Photron) with a recording rate ranging 
from 1500 to 30 000 fps (frames per second). The impact 
velocity and the droplet size were determined from the 
captured images. 

Several control parameters influence the impact dy- 
namics, for example, droplet size, liquid viscosity /i, im- 
pact velocity Vi. We describe these effects in terms of di- 
mensionless numbers: the Weber number We, the ratio 
of kinetic energy to surface energy, characterizing the de- 
formability of the droplet; the Reynolds number 7£e, the 
ratio of inertia to viscosity effect: 

We=^t, ne=^. (1) 
a /i 

Here R is the radius of the liquid drop, Vi is the impact 
velocity, p is the liquid density, a is the surface tension, 
and p is the liquid viscosity. In this paper, we did not 
simultaneously change We and 7Ze by using different liq- 
uids; instead, these parameters are only different ways 
to non-dimensionalize the velocity Vi of a milli-Q water 
droplet. The Ohnesorge number, comparing viscous and 
capillary forces, Oh = VWe/7^e = pj \J pRcr is small 
~ O(10 -3 ) in our experiments. In addition, we control 
and measure the pressure of the surrounding air, noted 
as Pair- For the data presented here we specify the cases 
under reduced air pressure, otherwise P a i r = 101.3 kPa 
at the standard ambient pressure. 

RESULTS AND DISCUSSION 

In this section, we compare impact events on CNFs 
and on microstructured surfaces with similar wettabil- 
ity. Fig. [2] shows a small fraction of the phase space 
and the time evolution of the impact dynamics for wa- 
ter drops impacting upon CNFJs of Osca = 155 ± 3°. 
In this small We-number regime (< 2.5), the impact 
events include non-bouncing Fakir droplet, complete re- 
bound, sticky wetting ball, and partial rebound, gener- 
ally in this sequence as We is gradually increased. When 
a water droplet is gently deposited, due to small kinetic 
energy the droplet maintains in a Cassie- wetting "Fakir" 
state with air trapped underneath the drop and thus a 
high contact angle during the whole impacting process, 
as shown in Fig. [2f. In this heterogeneous state, air is 
trapped underneath the droplet and between the porous 
fibers. The non- wetting nature of the sample limits the 
spreading of the drop and sometimes even leads to a com- 
plete rebound (Fig. Ek). At high Vi > 0.28 ms" 1 a wet- 
ting transition can occur as the kinetic energy overcomes 
the surface energy, associated with the liquid surface ten- 
sion, and the droplet turns to the completely wetting 
Wenzel state, in which it is pinned on the surface (Figs. Eh 
and d). 
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FIG. 2: (a) Fraction of the phase space of the impact dynamics upon a carbon nanofiber substrate (Fig. [TJi) with a static 
contact angle of (155 =b 3)° at small We numbers, with the time evolutions of the impact events shown in (c)-(f). (b) shows 
the corresponding impact velocity V% vs. the inverse contact time. 



The average contact time before the droplet bounces 
off the surface is 12.5 ms (based on the marked data 
(•) in Fig. (2}d). This average value is consistent with 
the value 14 ms found by a recent study with a multi- 
walled carbon nanotube array with a static contact angle 
of 163° [29j]. The coexistence of different dynamical be- 
haviors for 0.28 ms -1 < V\ < 0.35 ms -1 illustrates the 
complexity of the impact process, showing the interplay 
of several parameters and indicating that other factors 
can affect the phase space. For instance, from our results 
we note that the complete rebound regime happens at 
higher Vi between 0.5 ms -1 and 1.1 ms -1 for a more hy- 
drophobic CNF substrate with a larger Osca= (163d=4)°. 

Fig. [3] shows the phase diagram for We between 1 and 
10, corresponding to high 7£e, when a water droplet im- 
pacts upon a CNF surface of Bsc a = (152 ±3)° shown in 
Fig. [lb- In this We-range, as discussed above, the impact 
scenario can display a complete bouncing or a pinning of 
contact line resulting in wetting. Also, trapping an air 
bubble or jetting can happen due to the development of 
an air cavity. Fig. [4] shows the detailed evolutions of 
(a) trapping an air bubble and (b) emitting a fast jet. 
For our We-range between 2 and 8, on impact a sur- 
face capillary wave is excited and the droplet deforms 
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FIG. 3: (a) Impact dynamics upon a substrate of CNFJs, 
(shown in Fig. [TJ)) with a static contact angle of (152 ± 3)°. 
The insets show the snapshots of the impact events: sticky 
ball at the Wenzel state (■), partial rebound (a), trapping of 
an air bubble (♦), and jetting (*). 





FIG. 4: Time evolutions of the impact events of (a) trapping an air bubble at Bsc a = (163 =b 3)°, We = 5.3 and lZe= 571, and 
(b) jetting at Osca = (152 ± 3)°, We = 7.7 and IZe = 756. 



like a pyramid around 1.5 ms. The oscillation of the 
surface deformation often makes a toroidal droplet pro- 
ducing a cylinder-like cavity; see the snapshot at 3.2 ms 
in Fig. 2k and that at 2.2 ms in Fig. SJd. Then the fol- 
lowing dynamics of the cavity restoration determines the 
subsequent behaviors: in (a) at a lower Vi the surface 
wave restores the top of cavity sooner and closes it up 
producing a residual void; in (b) at a higher Vi the fast 
collapse of the cavity produces a fast jet. Our data reveal 
that the shooting jets can be as fine as 30 /im in radius 
and as fast as 9 x Vi. Similar observations underlying 
the same mechanism have recently been investigated with 
ordered, superhydrophobic microstructures for the We- 
number ranging between 0.6 and 16 [6[. The analogy of 
the impact phenomena suggests that the formation and 
the collapse of the air cavity may be insensitive to the 
details of surface roughness. 

In the higher We-number range (between 90 and 140) 
splashing impacts occur as shown by Fig. [5] for CNFs. 
Here, the "splashing" phenomenon refers to the forma- 
tion of many satellite droplets, which merge during the 
spreading and/or contracting stages of the water film 
from our observations. In Fig. [5^ under 1 atm, at 1 
ms tiny droplets radially emit while the water sheet still 
spreads outwards with wavy perturbations, resembling 
"fingers" , at the edge of water film. Several droplets orig- 
inated from the wavy perimeter develop between 1.8 ms 
and 5.8 ms while the main water sheet contracts, shaping 
into a partial rebound. A similar splashing phenomenon 
has lately been observed at We ^160 with a superhy- 



drophobic surface consisting of regular micro-pillars of 
h = 37 /im in height, and a = 3 /im in the interspac- 
ing distance between pillars 0- This Ref. also reports 
that the critical We-number above which droplets eject is 
about 700 for a flat hydrophobic solid and no fragmenta- 
tion of droplets even at We = 1000 for a flat hydrophilic 
solid. With these observations, the authors conjecture 
that the role of the film of air plays a crucial role because 
the non-wetting Fakir state promotes an air film for su- 
perhydrophobic solids. However, our investigations find 
the effect of air pressure on the splashing to be negligible 
in this parameter regime; Fig. [5b reveals the same kind 
of splashing at a low air pressure, P a { r = 14.7 /cPa, with 
similar We and IZe to those in Fig. [5^. Our experiments 
at different P a i r with similar We numbers resemble the 
same type of splashing. Hence, the data suggest that 
the mechanism of aforementioned splashing behaviors is 
due to the effect of surface roughness, and not due to the 
surrounding air. 

Fig. [6] shows a partial rebound for an impacting drop, 
with similar control parameters P a i r = 101.3 kPa, We = 
120 and IZe = 2882, but now for a polymeric microstruc- 
tured surface consisting of round pillars of a = 1.27 /im, 
w = 5 jim and h = 10 jim. As observed from the snap- 
shots, the rim is not destabilized in the same way as 
for the unstructured CNFJs (compare Fig. [5] and [6] at 
t ~ 1 ms). Instead of perimetric fingers, wavy bumps 
develop and eventually merge together at the center to 
form an elongating water fountain. The pinning of the 
water wets the microstructured surface and breaks the 
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FIG. 5: Snapshots of the splashing evolutions on CNFs of Osca — (163 =b 3)° at different air pressure (a) Pair — 101.3 kPa, 
We = 115.3 and lie = 2784, and (b) P air = 14.7 kPa, We = 141.7 and Tie = 3060. The inset bars indicate a length scale 
of 2 mm. 
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FIG. 6: Time evolutions of impact dynamics upon an ordered, microstructured polymeric substrate of Osca — (147 ±3)° at 
Pair — 101.3 kPa, We = 120 and He — 2882. This surface is comprised of round pillars of a = 1.27 /xm, w — 5 /xm and 
h = 10 fim arranged in a rectangular lattice and yields the roughness IZf = 5. 



fountain while surface tension is in action, resulting in 
a few breaking- up droplets. Consistent with the results 
in Ref. 0, our experiment at We = 147, He = 3524 
and P a i r = 101.3 kPa shows that a few (~ 5) satel- 
lite droplets develop during the contraction of the water 
sheet upon a micropatterned substrate of a = 5 //m, 
w = 5 //m and h = 20 //m with a roughness Tlf= 2.8. 
The comparisons between these different experiments re- 
veal that roughness has a profound effect on splashing 
impacts. In addition, the multiscale, uncontrolled sur- 



face roughness of CNFs display pronounced formations 
of satellite droplets, comparable with the impacts upon 
ordered micropatterned substrates under alike parame- 
ters. This finding suggests that the smaller length scale 
of carbon fibers may enhance the instability of perimetric 
wavy fingers and thus produce stronger fragmentations. 

Fig. displays the phase space of water impact dynam- 
ics onto polymeric microstructures of Osca — (147 ±3)° 
explored at different air pressure P a ir- Different col- 
ors depict different surface roughness IZf controlled by 
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FIG. 7: Phase space of water drop impact onto polymeric 
microstructures. Different colors indicate different roughness, 
IZf = 1.9 (green), 2.6 (magenta), 3.0 (black), and 5.0 (cyan), 
determined by the geometric arrangements of the micropillars 
of h = 10 /im. Different symbols present various dynami- 
cal impact behaviors: jetting (*), complete rebound (o), par- 
tial rebound (A), sticky ball (□), formation of a few satellite 
droplets while the main impact reveals a partial rebound (>) 
or a sticky wetting ball (<). 



CONCLUSIONS 

In comparison with the drop impacts performed upon 
microstructured hydrophobic surfaces, for We < 10 and 
under the atmospheric pressure impacts upon CNFs of 
multiscale roughness reveal similar dynamical behav- 
iors including Cassie- wetting Fakir droplet, complete re- 
bound, partial rebound, jetting, and trapping bubbles. 
This comparison implies that impact evolutions are in- 
sensitive to the details of the roughness of superhy- 
dropohobic surfaces at small We-number. Rather than 
a bouncing-pinning transition, driven by a threshold of 
impact velocity, as described in Ref [7[ for regular mi- 
crostructures, we noticed the coexistence of complete re- 
bound and sticky wetting droplet for CNFs at 1.1 < 
We < 1.9. This observation indicates that the details 
of roughness can influence the transitional boundary in 
the phase space of impact events. 

Partial rebounds occur in the We-range between 10 
and 130 for water drops hitting upon micropatterned sub- 
strates with h = 20 /im and a roughness ranging from 
2 to 10. Air pressure and the surface roughness have 
hardly any effect in this partially rebouncing regime. At 
higher We > 120, profound splashing impacts forming 
several satellite droplets take place for superhydrophobic 
CNFs. This splashing mechanism is not suppressed by 
the decrease in air pressure, which was found to debilitate 
the corona splashing on smooth wetting surfaces Io|, EH . 
Beyond We > 140, in contrast to the violent splashing 
events on the rough CNFJs, on imcrostruture substrates 
the fragmentation is much less pronounced, revealing the 
importance of sub-micron roughness for splashing. 
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